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ABSTRACT 

Observations of galaxy-galaxy lensing from Sloan Digital Sky Survey (SDSS) are com- 
bined with the TuUy-Fisher and fundamental plane relations to derive constraints on 
galactic halo profiles. We show that both for early and late type galaxies around 
the rotation velocity decreases significantly from its peak value at the optical radius 
to the virial radius r2oo, Wopt/^^200 ~ 1-8 with about 20% uncertainty. Such a decrease 
is expected in models in which the halo profile is very concentrated, so that it declines 
steeper than isothermal at large radii. This large decrease can be explained as a re- 
sult of both a concentrated dark matter profile and a significant stellar contribution 
to the rotation velocity at the optical radii. We model the stellar component with a 
thin rotationally supported disk or a Hernquist profile and use adiabatic dark matter 
response model to place limits on the halo concentration as a function of the stellar 
mass to light ratio. For reasonable values of the latter we find concentrations C200 con- 
sistent with CDM predictions, suggesting there is no evidence for low concentrations 
for the majority of halos in the universe. We also discuss the origin of Faber- Jackson 
relation L cx cr^ in light of L oc w|5q relation found for early type galaxies above 
from galaxy-galaxy lensing. This leads to a decrease in fopt/i'200 with luminosity 
above L*, so that at 7L* the ratio is 1.4. This is expected from the fundamental plane 
relation as a result of a reduction in the baryonic contribution to the total mass at 
the optical radius and a decrease in optical to virial rotation velocity in dark matter 
profile. These results imply that relations such as TuUy-Fisher and Faber-Jackson are 
not simply those between the mass of dark matter halo and galaxy luminosity, but are 
also significantly influenced by the baryonic effects on the rotation velocity at optical 
radii. 
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1 INTRODUCTION 

The study of dark matter profiles around galaxies has been 
an active area since the original discovery of dark matter 



from the flat rotation curves in spirals (see sofue & Rubin 
2001 for a review). These show that there must be dark 
matter in the outer parts of the galactic halos, but its ex- 
tent is uncertain because of the limited range probed by ob- 
servations. More recent observational studies of rotational 
velocity in dwarf and low surface brightness galaxies have 
suggested that the amount of dark matter in the central re- 



ulation bv CDM models (e.g. 


Moore 1994, 


Florcs & Primack 


1994, 


Debattista & Sellwood 200(\, 


do Blok ct al. 2001 


). The 



CDM models predict very concentrated dark matter haloes. 
This is usually parametrized using a universal dark matter 
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profile such as NEW ( jNavarro, Frenk, fc White 1997] ) f\. For 
a given virial mass of the halo, in this paper defined as the 
mass enclosed within a sphere of radius r2oo within which the 
density is 200 times critical density, the profile has a char- 
acteristic shape which depends only on a single parameter. 
For NFW profile, in which the slope continuously changes 
from the inner value of -1 to the outer value of -3, the free 
parameter is often defined as the concentration C200, which 
is the ratio between the radius r2oo and the scale radius 
Ts where the slope is close to —2. Higher concentration pa- 
rameters imply higher densities at the scale radius r^. On 



t Although other profiles have been proposed that differ signif- 
icantly from NFW in the inner par ts of the halo, the y agree 
well with NFW in the outer parts (Klypin ct al. 2001). Since 
for the present work we are concerned predominantly with the 
outer parts we only use NFW profile. 
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galactic scales CDM models predict C200 ~ 7 — 15 depending 
on the halo mass, matter density, shape and amplitude of 
the power spectrum. 

However, observational case for low concentrations is 
not cle ar and different conclusions have been reached by 



other studies (e.g. van den Bosch & Swaters 2001). Even 



if the observational results of low dark matter density in 
the inner regions are confirmed, they do not necessarily in- 
dicate a problem for the CDM models. The variation in 
the profile shapes between halos is large and it is not clear 
whether the dwarf or low surface brightness galaxies, which 
show strongest evidence for low density cores, can be as- 
sociated with the average halo population. For example, if 
these galaxies are associated with a population that had a 
major merger in the recent past, this may lead to a sig- 
nificantly flatter and less co ncentrated halo struc ture than 
the average population (e.g. Wechsler et al. 2001). Another 



possibility is that astrophysical processes such as bar rota- 
tion r edistribute the dark ma tter in the inner parts of the 
halo (Weinberg & Katz 2001). However, if the case of low 
concentrations is extended to the population of galaxies as 
a whole and to the outer parts of the halo, which cannot 
be affected by astrophysical processes, then the CDM crisis 
would become much more severe. 

It would thus be useful to have some information on the 
halo structure for the mean population of galaxies. While 
many observationally determined rotation velocity profiles 
exist, the main uncertainty has always been the relative 
contribution between the baryonic component in the gas, 
bulge/ellipsoid and disk and the dark matter in the halo. 
The baryonic and dark matter components are difficult to 
separate because the conversion from the star light (and, to 
a lesser extent, gas density) to baryonic mass is still uncer- 
tain from the stellar population synthesis models and other 
studies. The situation is further complicated by the fact that 
baryonic component is not dynamically negligible and dur- 
ing its condensation it induces a response of the dark mat- 
ter halo changing its original distribution in the inner parts. 
While these problems are ameliorated for the low surface 
brightness galaxies, they are almost never negligible. 

Since the modelling of relative disk or spheroid and halo 
contributions is rather difficult in the inner parts it is useful 
to con centrate on the outer parts of the halo, where the 



baryon influcTice is loss iinportant. Hnwnvnr. optical rntatigji. 



curves ivpically extend only out to 10 — 20kpc and even the 
ones based on HI measurements do not extend beyond 30- 
50kpc. Similarly, velocity dispersion studies in early type 
galaxies also do not extend past a few effective radii and 
even X-ray studies of large ellipticals do not ext end beyond 
10 effective radii (Loewenstein & White 1999). Over this 



limited range the data on disk L« galaxies indicate that 
the rotation curves rise in the inner parts of the disk and 
then stay approximately constant or decline somewhat out 
to the outer limit of observations. This flatness indicates that 
approximating the density profile as isothermal, p{r) oc r~^, 
is a good approximation to the matter profile over this range. 
On the other hand, theoretical CDM models predict that 
the velocity profile of the dark matter increases out to twice 
the scale radius and then slowly declines beyond that 
as the slope gradually decreases towards -3. The fact that 
we do not observe a decrease in rotation velocity at radii 
below 2.15rs ~ 30 — 50kpc is presumably due to the baryonic 



effects, which increase the mass in the center. The decrease 
at large radii is however a robust prediction of these models. 

While, until recently, no accurate determination of the 
mass profile at large radii has been available, r ecent galaxy- 



galax y lensing observations by the SDSS team (McKay et al 



2001) have improved the situation significantly by obtaining 
the morphology and luminosity dependence of the signal. 
Theoretical analysis of lensing data must take into account 
not just galactic halos, but also those from groups and clus- 
t ers, which dominate on large scales above 200-300/i~^kpc 
( [Cuzik fc Seljak 2002 ). These results show that a late type 
L, galaxy, with Li — 2.7 x 10^° Lq (where we have ap- 
plied a 30% internal extinction correction in I band), has 
a mass M200 ~ (3.4 ± 2.1) X lO"/i"^M0, corresponding 
to the circular velocity at the virial radius U200 '~ 115km/s. 
While the error is still quite large it is gaussian distributed in 
mass, so an increase in W200 by 35% to 150km/s is excluded 
at 95% confidence level. This number can be compared to 
the maximum rotation velocity for such a galaxy, wh ich is 
208km/s with a small scatter (Giovanelli et al. 1997), the 
well known Tully-Fisher relation. Comparing the two values 
shows that a decrease in rotation velocity from the optical 
to the virial radius indeed occurs for the average popula- 
tion of spiral galaxies. The decrease is large, almost a factor 
of 2 and even if we push the virial velocity up by 2-a to 
^200 = 150km/s the decrease is still around 40%. As shown 
in this paper, a similarly large decrease in rotation veloc- 
ity is obtained also for early type galaxies which are not 
rotationally supported. 

These results are interesting, since they are in the di- 
rection predicted by the CDM models and demonstrate that 
the halo profiles indeed become steeper than -2 in the outer 
parts of the halo. They set tight limits both on the structure 
formation models, by limiting the acceptable range of con- 
centration parameters, and on the disk/spheroid formation 
models, by constraining the stellar mass to light ratio. Only 
for specific values of stellar mass to light ratio and halo con- 
centration can one satisfy these constraints. The purpose of 
this paper is to investigate the constraints in detail using 
halo and disk/spheroid formation models. 

Previous work on this subject has explored the con- 
straints from the rotation velocity at optical radius given 
by the zero point of Tully-Fishe r relation ([Mo 



Mao 



200C 



Navarro fc Steinmetz 2000| , |Eke, Navarro, fc Stein- 



metz 2001). In the absence of virial mass information hese 



models must rely on additional assumptions to derive the 
constraints on cosmological models. The advantage of the 
additional information from lensing is that it provides an- 
other dynamical constraint at large radii, which can remove 
some of the modelling uncertainties present in previous mod- 
elling. In addition, while previous work only explored the 
constraints from late type galaxies, in this paper we also 
investigate the constraints from the early type galaxies. We 
find these are more robust both because the virial masses are 
more accurately determined and because the velocity dis- 
persions at optical radii are obtained from the same SDSS 
sample. 
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2 LATE TYPE GALAXIES 

The average rotation velocity at optical radii can be ob- 
tained from the Tully-Fisher (TF) relation, which in /-band 
is given by (GiovancUi et al. 1997) 



Li = 2.7 X 10' 



(208km/s 



Vopt 



-0' 



where we used I — Slog/i = —21 — 7.68[log(2«, 



'opt J 



(1) 

2.5] 



and 7, 







4.15. We denote with Vopt the maximum rotation 
velocity typically achieved at the optical radius (roughly 3 
times the scale radius of the disk Rd, which is not to be 
confused by the scale radius of the halo r^). At this radius 
the rotation curve still has a significant contribution from 
the disk. Rotation curves at larger radii show that, for this 
range of luminosities, the rotation curve at L« are flat or 



decline blightly out t o the largest radius observable, typically 



a few o ptical radn (^scrtano & van Gorkom liJiJl Persic 
Salucci j fc Stel IQQq [Verheijen 200l'^ 

We would like to compare the dark matter velocity in 
the inner parts of the halo to the SDSS galaxy-galaxy lensing 
results. For a late type L, galaxy with i* — 5 log h = —21.26 
one finds M200 = 3.4 x 10"fc-^ Mfi^(c20o/10)~" " (iGuzik & 
Seljak j002[ ). Using Iq = 4.52 ( [Blanton et al. 200]] ), apply- 
i ng an average .3 magnitude internal extinction correction 
(Verheijen 2001) and converting from the viriai mass to the 



viriai velocity using the relation GM200 /t200 = 1^200 ^6 find 



: 2.7 X 10' 



i 115km/s 



«200 



/ C2OO 

vw 



0- 



(2) 



The dynamical range of galaxy-galaxy lensing is still rather 
small and for late type galaxies luminosity dependence of the 
viriai mass cannot be established using the present sample, 
so we focus on L* galaxies, which dominate the late galaxy 
galaxy-galaxy (g-g) lensing signal. 

By combining the two equations above we can deter- 
mine the best fitted ratio 

fopt 

W200 

Since the viriai mass of late type galaxies is consistent with 
at 2-a level there is no upper limit to the velocity ratio, while 
the lower limit is given by Dopt/u200 > 1.4 at 2-(t, where the 



1.8. 



(3) 



error bi|idgeL includes lli u slalislical uiiuis on l liu zuiu puiiii 
of TuU) 



Fiijhtii" relation (uiovancUi et al. 1007^7l)ut m dom 



inated by Llie eiiui un Llie v iiial mass fiuni galaxy-galaxy 
lensing (Guzik & Seljak 2002) and we ignored the small con- 



centration dependence. In addition to the statistical errors 
there are also possible systematic difi'erences b etween the 



different T ully-Fisher zero po i nt determina tions ( Giovanelh 
et al. ld97|, IWillick et al. 19961, fehanks 19971) which can be up 



to 0.2 magnitude. Similarly there could be svstem atic dif- 
ference s be tween color selected late type galaxies ( McKay 
et al. 2 )01) and those selected for rotation velocity studies, 
although morphological studies indicate that the late type 
sample is domi nated by Sb/Sc mor phological type used also 
in TF studies ( |Strateva et al. 200l| ). 

What constraints does equation ^ imply on the struc- 
ture formation models? As discussed above, a decrease in 
rotation velocity implies that the mass profile is steeper 
than isothermal in the outer parts of the halo. This can 
be either because the dark matter profiles are steep or be- 
cause the stellar disk has a significant contribution to the 



rotation velocity (or both). An NFW profile is given by 
p{x) = psX~^{l + x)~^ , where x = r/r^. The rotation ve- 
locity is VcifY = GM{r)/r, where M(r) = 47vpsrgA{c) 
and A(c) = ln(l -I- c) - c/(l -I- c). For NFW profile the 
rotation velocity increases at small x up to the peak at 
X = 2.16 and then declines gradually to the viriai ra- 
dius. The ratio between the maximum and viria i velocity 
is Umax/u200 = 0.46 [c/y4(c)] jBuUock et al. 200l| ). For the 
range of concentration parameters predicted by CDM mod- 
els (c < 20) one finds Wopt/«200 < 1-4. It is thus unlikely that 
the dark matter halo can explain the decrease by itself, since 
even for very concentrated halos the velocity decrease is less 
than observed. We must therefore include the contribution 
from the disk. 

We will use a model with a thin exponential disk in a 



dark matter halo ( Mo, Mao, fc White 1998 



van den Bosch 



200C). We model the disk as a thin exponential surface den- 



sity profile, 

E(i?) = Eoexp{~R/Rd 



(4) 



with the total disk mass given by Md = 27rEo-Rd- The disk 
mass is related to the disk luminosity using disk mass to 
light ratio Tj — M/L. We will neglect the bulge contri- 
buti on both to the luminosity and to the rotation curve 
(see Mo et al. 199^ for a discussion of this assumption). 
We will use observations to provide the typical scale length 
Rd of the galactic disks. For an L, galaxy obse rvations give 
Rd ~ 3.5/i"^kpc ( ICourteau 1997| , |de Jong 199^ ). The scatter 
for this quantity at a given luminosity is rather large, since 
galaxies come with a range of surface brightnesses. Theoret- 
ically, the scatter ha s been linked to the scatter in the spin 
parameter of haloes (Mo et al. 1998). However, since in this 



paper we are primarily concerned with the average proper- 
ties of galaxies and not in the scatter around the mean we 
will assume Rd — 3.5ft~^kpc in the analysis as a mean value 
for a typical L« galaxy. 

Disk gravity contributes to the measured rotation ve- 
locity. In addition, disk gravity also induces a response of 
the dark matter halo in the inner regions. The standard ap- 
proach to model this is to assume an adiabatic contraction 
of the halo, which remains spherical, so that the angular 
momentum of individual particles is conserved (see 



White 1984 



Blumenthal et al. 1986 and 



Barnes 



Navarro & Stein 



mctz 200C for more details and numerical tests on the valid- 



ity of this model) . This leads to an implicit equation for the 
final radius rf of dark matter mass as a function of initial 
radius ri 



AfDM(n)n = [MDM(n)(l - /,) + Me(rf)]r-f, 



(5) 



where Mdm('") and Ms{r) are dark matter and stellar mass, 
respectively and /* is the stellar mass fraction of the halo. 
Assuming that the baryons which do not end up in the disk 
have the same distribution as the dark matter allows one to 
solve the system completely for a given disk mass and scale 
length and for a given halo profile. Flattened nature of the 
disk is used when obtaining stellar rota tion velocity from the 
mass profile (Binney & Tremaine 1987). As a free parameter 
we will use the stellar mass to light ratio T/ (expressed in 
solar units), which from the known viriai mass and lumi- 
nosity of Lt galaxy (and ignoring bulge contribution to the 
luminosity) can be related to /* as 
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Ti = Md/L, = 



3.4 X 10"/i"^Mq/* 



^0 



12/* 











(6) 



= (1 - 2)h (e.j 

Note that /, should not exceed 



Bottema 1997) 



Typical values are Tj 
giving /. = 0.1 - 0.2 
Qb/^m ~ 0.04/f2m, since only the baryons within the virial 
radius can condense to make stars. This implies Qm < 0.4 
for the fiducial value of M, and Qm < 1 for the 95 % c.l. 
on M«. This is in agreement with other determinations that 
give Qm < 0.4, so the fraction of baryons converted to stars 
does not exceed the available supply, although it comes quite 
close to this limit and is an argument against high values for 
Tj. 

We can solve for Hopt 7^^200 for any given /, (or T/) and 
C200. An example with T/ = 1.7h and C200 ~ 12 is shown 
in figure ^ One can see that the velocity ratio of 1.8 can be 
naturally obtained in a model with concentration pararneter 
in the r ange predicted by CDM models, 8 < 0200 < 15 (Eke 



Navarro, & Steinmetz 2001, Bullock et al. 2001), and with 
the expected stellar mass to light ratio, Ih < T i < 2h. The 
rotation velocity is reasonably close to flat over the optical 
region, but decreases by 10-20% out to the largest range 
observable in HI, in agree ment with the observations for 



this range of lu minosities (Casertano & van Gorkom 1991 



Verheijen 2001). The disk and dark matter contributions 
to the mass within the optical radius are comparable, each 
contributing 50% in this example. This implies that this 
model satisfles the requirement that the zero point of TF 
relation is independent of the disk surface brightness, which 
requires about 50% dark m atter contribution to the rotation 
velocity at optical radius (Mo & Mao 2000). Note that the 



adiabatic response of dark matter is quite signiflcant and 
dark matter would be subdominant if it were not compressed 
by baryonic condensation. Because of this the dark matter 
contribution to the rotation curve is never negligible, even in 
the inner parts of the galaxy. However, our predictions may 
not be reliable inside the optical radius, where bulge makes 
a significant contribution and adibat ic approximation may 
not be vaUd ([Weinberg fc Katz 200"l[). 



The example above was chosen based on the typical val- 
ues for concentration and stellar mass to light ratio. More 
generally for any choice of one there will be a particular value 
for the other that satisfies the constraint in equation 
This is shown in figure |^ for a family of stellar mass to light 
ratios. Low values of T/ cannot be made compatible with 
the constraint in equation ^ unless the concentrations are 
unreasonably high (e.g. c > 20 for T/ < Ih). Low concentra- 
tions are also not acceptable unless the observed i'opt/w200 
is decreased by 2-a to 1.4. In this case one must adopt ei- 
ther very low concentrations for reasonable T/ ~ 1.5h or 
very low T/ < 0.7h for reasonable c ~ 10. On the other 
hand, a positive deviation in i'opt/i'200 can be explained by 
using a somewhat higher values for T/. We conclude that 
the results are just what is expected from the CDM models 
with standard concentrations and standard stellar mass to 
light ratios. Only if the virial masses deviate by 2-a in the 
positive direction from the mean value does one run into 
problems with the standard stellar mass to light ratios and 
one requires T/ = 0.5h. The majority of late type galaxies 
therefore show no evidence for shallow density profiles in the 
outer parts of the halo. 
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Figure 1. Rotation velocity contributions from the disk (short 
dashed), dark matter (dotted) and total (solid) out to the virial 
radius. Also shown is the dark matter without adiabatic response 
to baryon contraction (long dashed). Insert plots the same over 
the radii accessible with HI data, showing that the rotation curve 
is close to flat over this range. We do not model the bulge/bar, 
so the predictions within inner few kpc are not reliable. 
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^ T\=2.5hJ 
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______ — " T,= 1.5hj; 
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Figure 2. Lines of constant T/ (times K) as a function of c ver- 
sus "Oopt /''200 . Thick dashed horizontal line is the observationally 
determined value fopt/f200 = 1-8- A. 2-u (98% confidence level) 
lower limit is Dopt/f200 > 1.4 (thin dashed line), while currently 
the data provide no upper limit since the virial masses are con- 
sistent with 0. 
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3 EARLY TYPE GALAXIES 

Early type galaxies have some advantages in constraining 
the outer parts of the halo profiles. The main advantage is 
that they show a stronger g-g lensing signal, so that both 
Mt and its scaling with luminosity are reliably determined 
from the current data and the errors associated with it are 
significantly smaller. Moreover, one can study their early 
type dynamics using the same spectroscopic SDSS sample 



also us ;d tor g-g lensing (the actual samples used in the 
analysis here are not complete ly equal since the an alysis 



for the early type dynamics in Bernardi et al. (2001) was 
using a lar ger sample than the lensing analysis in McKay 
et al. (2001 ), but the statistical properties of the two samples 



should be ver y similar). The fundam ental plane relations 



as derived by Bernardi et al. (2001) are L oc a (x R, 



where a is the central stellar velocity dispersion and _Rc is 
the effective radius of de Vauculeurs profile. The values at 
L. in i' = -21.26 are a = 177km/s and = 2.7/i"^kpc 
(note that we are usin g L, for the luminos ity function of the 
whole sample given in Blanton et al. 2001, not just the early 
type). 

For early type galaxies the conversion from the stel- 
lar velocity dispersion to matter circular velocity is less 
straightforward. In principle one can obtain it by solving 
the Jeans equation, which however depends on the unknown 
anisotropy of velocity dispersion. If the circular velocity of 
the matter does not change very much over the optical re- 
gion then virial theorem guarantees that iiopt ~ 3^^^ (a), 
if (a) is luminosity weighted line of sight velocity disper- 
sion. In practice luminosity average is difficult to achieve 
and the velocity dispersion decreases with radius, so more 
often observers report the central velocity dispersion (Tcontrai ■ 
Based on detailed kinematic analysis and on comparison 



betwe en strong lensing and stellar dispersions Kochanek 
(1994) argues that in this case the relation is closer to 
Vopt = 2^^'^acentrai. In SDSS the aperture is determined by 
SDSS fibers and is rflbor = 1.5", implying that for nearby 
galaxies only central parts of the galaxy are detect ed, while 
for distant g alaxies most of the light is observed. Bernardi 



et al. (2001 ) attempt to correct for this using an empirical 
fit, CTcentrai/o" = (Srflbor /^c , to obtaiu the central veloc- 
ity dispersion. The correction is empirical and appears to be 
somewhat small to account for the suggested 22% difference 
between the central and luminosity weighted velocity disper- 
sion, since one needs very small i?c = 0.15" to achieve this. 
We will use an intermediate value of Vopt ~ 1.5a to convert 
from the central velocity dispersion to the rotation velocity 
at the optical radius. This should have at most a 10% sys- 
tematic uncertainty attached to it. This conversion agrees 
well with the studies of slowly rotating elliptical galaxies 
where both Vopt and cTccntrai have been measured (Gerhard 



et al. 200l[ ). 

Another ingredient in the modelling is the dark matter 
response to baryonic contraction. We model it again using 
the adiabatic model. In principle there is no reason why 
such a model would be appropriate for elliptical galaxies, 
where stars are not on circular orbits, but numerical simu- 
lations of galaxy formation have found that adiabatic com- 
pression model works remarka bly well even for such systems 
(Gotbrath & Steinmetz 2001), perhaps as a consequence of 



Hernquist profile ( iHernquist 199C| ), which has been shown 
to give a light profile very close to the de Vaucouleurs profile 
and has an a nalytic 3-d rad ial distribution, to model the star 
distribution ( Keeton 2001 ). Our canonic value for the stel- 
lar mass to light ratio in i' is a factor of 2 higher than that 
of late type galaxies in the same band, T^/ = SHMq/Lq. 
This is based on K — i' color difference between early and 
late type galaxies in SDSS, which is around 0. 2-0.3 magni 



tudc s (without internal extinction correction, Ivezic et al 



2001) and the fact that in K band luminosity to stellar mass 



conversion only depends on the age of population and dif- 
fers by less than a factor of 2 between early and late type 
galaxies for reason able IMF and assuming ages above 3Gyr 
(Drory et al. 2002). This stellar mass to light ratio is again 
rather uncertain and significantly higher values have been 
suggested in the literature. Most of the direct studies are 
done in B band, but even after correcting for a factor of 2-3 
difference between B and i' luminosity for early type galaxies 
it is on the low side based on the dynam ical studies of centra l 
regions using the minimal halo models (Gerhard et al. 2001). 



The simplest explanation is that part of the mass is actually 
due to dark matter, as discussed further below. Note that 
for Til = 4hMQ /Lq the stellar to virial halo mass ratio for 
early and late type galaxies become equal and are approach- 
ing the maximal baryon to dark matter ratio still allowed by 
the observations. It is thus unlikely that the average stellar 
mass to light ratio can significantly exceed this value if the 
virial masses from g-g lensing are correct. 

At L, = 2 X 10^"/i~^Lq the vir ial mass for early typ e 
galaxy is (9.3 ± 2.2) x lO"/i"^M0 ( |Guzik fc Seljak 2002D , 
which translates to V2oo ~ (160 ± 15) km/s at r2oo ~ 
160/i~^kpc. At the effective radius Rc — 2.7/i~^kpc the 
rotation velocity from optical velocity dispe rsion is fopt ''^ 
1.5 X 177km/s = 265km/s with a small error ( Bernardi et al 



2001), leading to 



Vopt 
V200 



;±o.2. 



(7) 



conservation of radial action in such systems. We use the 



The error is dominated by stellar to dark matter velocity dis- 
persion conversion and virial mass uncertainty. Here again 
we have a very large decrease from the optical to the virial 
radius, which is inconsistent with the flat rotation curve at 
more than 3-a level. This decrease is similarly large to the 
one observed for the late type galaxies. Such a decrease can- 
not be explained by the dark matter alone unless halos are 
extremely concentrated. More realistic models must include 
baryons, which make a significant contribution to the ro- 
tation velocity at optical radii, both by direct contribution 
and by compressing the dark matter. 

The resulting velocity profiles are shown in top of fig- 
ure ^ for the canonic values T/ = 3/i and C200 = 10 (we use 
a somewhat lower C200 than for late types since concentra- 
tion is expected to decrease with halo mass). The maximum 
rotation velocity peaks very close to the optical radius and 
has a value of Vmax = 270km/s. This is in a close agreement 
with the observed value Vopt ~ 265km/s and is well within 
the estimated error of 30km/s, indicating that this model 
has no problem explaining the observed ratio of the optical 
to virial rotation velocity. Note that at the optical radius 
the baryon and dark matter contributions are comparable, 
while at somewhat larger radii dark matter dominates. The 
resulting profile is much closer to a constant velocity SIS 
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Figure 3. Rotation velocity contributions from the spheroid 
(short dashed), dark matter (dotted) and total (solid) for an L» 
early type galaxy. Also shown is the dark matter without adia- 
batic response to baryon contraction (long dashed). Insert plots 
the same over the radii accessible with optical data, showing that 
the profile is close to isothermal over this range. 



profile than if just ligiit was contributing to tlie mass. Tfiis 
is in a good agreement witfi tfie conclusions from strong 
lensing ( Kochanek 199?:), r otation velocity studies of ellip- 
t icals ( Gerhard et al. 2001) and X-ray studies of ellipticals 
(Loewenstein & White 1999), which indicate that flat rota- 



tion curve over the optical region (r < 10/i~^kpc) is a better 
fit to the data than a constant stellar mass to light ratio with 
no dark matter contribution. 

A general exploration of c versus Ti is shown in figure 
^ for L*. Within 1-cr of the best fitted value for i'opt/i'200 
one has 2h < Ti < 4/i assuming c = 10, whereas at 2-cr level 
this is extended to Ih < Ti < 5h. While the range of stel- 
lar mass to light ratios suggested in the literature is rather 
large and extends even above Ti > 5h, such high values 
are typically found for minimum halo models and are thus 
an upper limit. Our analysis suggests that Ti < 5h both 
because of the dynamical constraint and because of the lim- 
ited baryon supply. Very high values of stellar mass to light 
ratios are also not compatible with the stellar population 
synthe sis model s using the observed IMF (e.g. 
& Chaillot 199^). 



Kauffmann 



S ince the scaling of vir ial velocity with luminosity L oc 



?;2oo ( G!uzik fc Seljak 2002 ) differs from the scaling of op- 



tical velocity with luminosity L oc Vopt above L* we must 
also compare at a higher luminosity. We choose L = 71/*, 
which corresponds t o the highest lumin osity bin in SDSS 



g-g lensing analysis (McKay et al. 2001) and is dominated 



by earl y type galaxies. At this luminosity one finds M; 



10"fc~ tM(7) (Puzik fc Seljak 2002D, which corresponds to 



^^200 = (3 10 ± 30)km/s. The corresponding central velocity 
dispersion is Ucentrai = 290km /s, implying «opt — 435km/s 
and Re = 10/i~^kpc (Bornardi et al. 2001). Here the ratio 



Figure 4. Same as figure Q for early type galaxies at L*. Also 
shown (thick dashed) is the best fitted value Vopt 7^^200 = 1-68 
together with 2-a (95% confidence level) lower and upper limits 
(thin dashed). 



is Vmax/v200 ~ 1-4 ± 0.2. The dominating error is the virial 
mass at this luminosity. 

Results of the adiabatic model calculation are shown in 
figure Is] using Ti/ = 3/i. We have used an even lower concen- 
tration value c — 8 for this hi gher mass halo, sinc e numerical 



simulations find c oc M'" '** ( Bullock et al. 2001 ). We agam 



find the rotation velocity at the optical radius exceeding that 
at the virial radius, but the excess is significantly smaller 
now. The rotation curve is very fiat for r > 5/i~^kpc with the 
value at the optical radius around 445km/s, in good agree- 
ment with the observed value of 435km/s. The peak value 
of the dark matter rotation curve alone is around 400km/s 
and is close to the value at the optical radius given that the 
rotation curve is so flat. Note that at this luminosity the 
dark matter dominates already at the optical radius and it 
is only below 5/i~^kpc that baryon mass exceeds that of the 
dark matter. However, the baryons do have a signiflcant ef- 
fect on the dark matter through the adiabatic compression, 
so that even though the dark matter dominates at the opti- 
cal radius it would have been comparable to baryons if there 
was no dark matter compression. Our analysis suggests that 
the data are consistent with brighter ellipticals being more 
dark matter dominated at the optical radii than the fainter 
ones and the stellar mass to light ratio not varying with 
luminosity. 

It is worth noting that the model can be extended by 
yet another order of magnitude in virial mass, to cluster 
masses with an elliptical cD galaxy at the center. A typ- 
ical example is Virgo cluster with M87 at the center. For 
M87/Virgo the rotation velocity increases fro m the optical 
radius te ~ 5/i~^k pc value Vopt ~ 450km/s (Romanowsky 



Kochanek 2001), to around 950km/s deduced from the 



observed X-ray temperature around 3keV. An example that 
can reproduce these constraints is shown in figure ^ where a 
halo with M200 = 2 x 10^^/i~'^Mq, c = 6 and spheroid with 
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Figure 5. Same as figure ^ for 7L* early type galaxy. 



Msteiiar = 2 X IO^/i^^A'/q have been used ( |Giraud 1999| ). 
In this case one has a ratio in rotation velocity from optical 
to virial radius i'opt/^200 ~ 0.5. This is caused primarily by 
the cluster halo profile becoming shallower than r~ in the 
inner parts. While baryons in cD increase the rotation ve- 
locity from the pure dark matter value in the inner parts of 
the halo they cannot pull in outer parts of the dark matter 
halo, which remain unaffected. The trend of optical to virial 
rotation velocity ratio decreasing with halo mass observed 
in elliptical galaxies is therefore expected as a consequence 
of more massive halos being less concentrated and less steep 
than isothermal in the central region, although for quanti- 
tative predictions baryons and their effect on dark matter 
must also be included. 



4 DISCUSSION 



Figure 6. Same as figure y for 2 X W^^h ^Mq early type galaxy 
in a 2 X lO^**h-^M0 halo as a model for M87/Virgo. 



row range of scales observed. This demostrates the power of 
g-g lensing measuring the mass at large radii combined with 
the more traditional methods that measure mass at smaller 

for a similar conclusion for 



radii (see also Wilson et al. 2001 



early type galaxies at higher redshifts). The concentration 
parameters obtained using reasonable stellar mass to light 
ratios are around c ~ 10 and are in a good agreement with 
predictions of ACDM model with i},n = 0. 3 and as = 0-8 , 
which gives C200 ~ 10 at M ~ 10^^ ft" ^Mq ( [Eke et al. 20011 ). 
We find no evidence for low halo concentrations in the main 
galaxy population both for early and late type galaxies. 

The theoretically predicted stellar mass to light ratios 
have considerable theoretical uncertainties. Assuming CDM 
halo profiles and concentrations we find T/ 1.5h for late 
type galaxies and Ti ~ 3h for early type galaxies. These 
values are in a good a greement with the stellar population 



In this 3aper we compare average properties of rotation ve- et al. 2002). While higher stellar mass to light rati os are 



synthesis m odels (e.g. Kaufi'mann fc Chariot 199q , Drory 



locities between optical and virial radii based on optical and 
galaxy-galay lensing measurements, respectively. Our model 
is statistical in nature, since it is not based on analysis of 
individual galajcies, but on their average properties as a func- 
tion of luminosity. On the other hand, the galaxies in our 
sample were not chosen on the basis of any selection crite- 
ria, so our results should apply to the galaxy population as 
a whole. Moreover, the large dynamical range between op- 
tical region (of order a few kpc) and virial region (of order 
a few hundred kpc) allows one to measure slow changes in 
rotation velocity which are not possible to detect from each 
of the observations individually. 

Our main conclusion is that the rotation velocity in 
galaxies decreases significantly from optical radii to virial 
radii. Such a decrease is theoretically expected since the dark 
matter profiles in CDM models are steeper than isothermal 
at large radii and contribution from stars further increases 
the rotation velocity at the optical radius. It has however 
not been clearly demonstrated previously because of a nar- 



often quoted in the literature (e.g. Serhard et al. 2001) these 



are based on minimal halo models and are thus often an 
upper limit. It is interesting to note that in our models the 
dark matter contribution to rotation velocity is comparable 
to that of stars even below efi'ective radius. This is because 
of the adiabatic response, which compresses the dark matter 
in the center. While this model may be oversimplified it 
nevertheless suggests that it may be difficult or impossible 
to separate the two components on the basis of dynamical 
studies. Using the minimal halo assumption the stellar mass 
to light ratio in early type galaxies may be overestimated 
up to a factor of two. 

While the optical to virial velocity ratio is above unity 
both for early and late type galaxies, it is decreasing with 
halo mass from 1.8 at 3 x lO"/i"^M0 to 1.4 at 10"fe"^Mr7 



^0 



(this statement is valid only for galaxies above L* and it is 
possible that the trend is reversed towards the low luminos- 
ity galaxies). This trend continues further into the cluster 
halo masses, where the ratio falls below unity. Such a trend 
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is expected in models where halo profiles are less concen- 
trated for higher halo masses, implying that the turnaround 
from an increase to a decrease in rotation velocity occurs at a 
larger radius relative to the virial radius. In addition, in more 
massive halos stars play a less important role both as a direct 
contribution to the rotation velocity and through their effect 
on the dark matter. The rotation velocity-luminosity scal- 
ings at optical radii, such as TuUy-Fisher and Faber- Jackson 
relations, are not directly related to the properties of dark 
matter, but also require a proper modelling of baryons and 
dark matter respons e to baryonic contraction (see also a 



relat ed discussion in [Gonzalez et al 



2000 



Kochanek 2001 



and van den Bosch 2001). While there are still uncertain- 
ties in the modelling of these processes, the simple models 
presented here reproduce well the constraints from the data 
both for early and late type galaxies. 

How do our results compar e to previous wo r k? The zero 



point of TP relation problem ( |Eke et al. 200l| , |Mo fc Mao 



200C)if the closest to the TF analysis done here. In the ab- 
sence of virial mass information the value of rotation velocity 
at a given luminosity does not suffice to make any general 
conclusions, so in general one has to make additional as- 
sumptions and/or modelling. For example, the stellar mass 
fraction in the halo /, obtained in previous work was lower, 
which lead to a higher virial mass for a given luminosity, 
which in turn requires lower concentrations and/or stellar 
mass to light ratios. By increasing /* close to its maximum 
value the virial mass can be reduced and this alleviates the 
problem. The same solution also solves the suggested over- 
predicti on of dark matter at the solar radius in our own 
galaxy (Eke, Navarro, & Steinmetz 2001), since again if the 



stellar fraction is higher the virial mass can be lower (there 
may be additional prob lems for CDM profiles in the inner 



parts of our galaxy; e.g. Binney & Evans 2001) 



For early type galaxies it has been suggested tha t the 
concentrations are low from the strong lens statistics ( Kee- 
ton 20c(l|), since very concentrated halos would overpredict 
the expected number of lenses. This is a difficult method 
to use since the expected number of lenses is very sensitive 
to the assumed luminosity function for early type galaxies 
as a function of redshift, which still has considerable uncer- 
tainty. Additional uncertainty arises from the adopted values 
for stellar mass to light ratios, which again can change the 
lensing statistics significantly. The lensing results are still 
compatible with low density CDM models, suggesting there 
is no discrepancy with our results, although more work is 
required to study this in detail. 

There are other problems that have been suggested as 
troublesome for CDM, such as detailed shapes of velocity 
profil es in the optical region (de B lok et al. 2001), bar rota- 
tion ( Debattista fc Se llwood 2000|) or t he halostructure of 
the Milky Way ( Binney fc Evans 2001 ). These probe inner 
regions of the galaxy where complicated physical processes 
may be taking place, so there is considerable more uncer- 
tainty in their theoretical predictions. For example, there are 
proces ses such as bar rotation that can disrupt dark matter 
cusps (Weinberg & Katz 2001). It has recently been shown 
that CDM profil es can fit most of the rot ation curves for 



normal galaxies (Jimenez, Verde, & Oh 2002). The galaxies 
that appear to be a problem for CDM belong to one of the 
specific subsamples, such as low surface brightness, dwarf 
or barred galaxies. It is possible, although not necessarily 



easy to arrange, that these samples are qualitatively differ- 
ent from the main population, for example by forming later 
and thus being less concentrated. Yet another possibility is 
that problems arise only below L* , since our analysis is only 
valid for galaxies around and above L*. Clearly, more work 
is required to resolve these issues. However, if the g-g lensing 
masses are correct, then for the main population of galaxies 
around and above L« the CDM model predictions for the 
amount of dark matter outside the inner few kpc do not 
exceed the observations, suggesting that the problems for 
CDM may not be as fundamental as previously suggested. 
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